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Abstract 

^ ■ We study the stability of /(i?)-AdS (Schwarzschild-AdS) black hole obtained from 

f(R) gravity. In order to resolve the difficulty of solving fourth order linearized equa- 
tions, we transform f(R) gravity into the scalar-tensor theory by introducing two 
\Q , auxiliary scalars. In this case, the linearized curvature scalar becomes a dynamical 

scalaron, showing that all linearized equations are second order. Using the positivity 

O ■ 

q-i of gravitational potentials and S-deformed technique allows us to guarantee the sta- 

^ny bility of /(-R)-AdS black hole if the scalaron mass squared satisfies the Breitenlohner- 

Freedman bound. This is confirmed by computing quasinormal frequencies of the 
CN , scalaron for /(i?)-AdS black hole. 
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1 Introduction 



f(R) gravities [TJ EJ [3] have much attention as one of strong candidates for explaining 
the current accelerating universe [I]. f{R) gravities can be considered as Einstein gravity 
with an additional scalar. For example, it was shown that the metric-/(i?) gravity is 
equivalent to the wbd = Brans-Dicke (BD) theory with a certain potential [S]. However, 
in order that f(R) gravities are acceptable, they obey minimal requirements for theoretical 
viability [21 |5]. Three important requirements are included: (i) they possess the correct 
cosmological dynamics, (ii) they are free from instabilities of tachyon and ghost |U [S] , 
(iii) they attain correct Newtonian and post-Newtonian limits. 

The Schwarzschild-de Sitter black hole was first obtained for a positively constant cur- 
vature scalar in f(R) gravities [9] and other black hole solution was recently found for a 
nonconstant curvature scalar in f(R) gravities [TU]. A black hole solution was also obtained 
from f(R) gravities by requiring the negative constant curvature scalar R = R < [TTj . 
For 1 + f'(R) > 0, this is similar to the Schwarzschild-AdS (SAdS) black hole obtained 
from Einstein gravity. In order to obtain the constant curvature black hole from u f(R) 
gravity coupled to the matter", the trace of its stress-energy tensor should be zero. 
Two solutions were found when the Maxwell [TTJ and Yang-Mills fields [12] are coupled to 
f(R) gravities. 

Gravitational stability of a black hole is a main issue of testing the adequateness of the 
solution [TSJ[T3]. However, the stability analysis seems not directly applicable to f(R) black 
holes including the Kerr black hole because f(R) gravity contains fourth-order derivative 
terms in the linearized equations [T5| HE] . In this case, the requirement (ii)[they are free 
from instabilities of tachyon and ghost] would come into play in testing the stability of f(R) 
black holes [17] . One may transform f(R) gravity into the scalar-tensor theory to eliminate 
fourth-order derivative terms by introducing two scalar fields [TS]. Then, the linearized 
curvature scalar became a scalaron, guaranteeing that all linearized equations are second 
order. For f(R) black hole, its linearized equations have led to the same equations for the 
BD theory [191 12H] • Using the stability analysis of black hole in the massive BD theory [20] , 
it was shown that the f(R) black hole is stable against the external perturbations if the 
scalaron does not have a tachyonic mass in asymptotically flat spacetimes [2~T] . 

In this work, we investigate the stability of f(R)- AdS (SAdS) black hole. We transform 
f(R) gravity into the scalar-tensor theory to avoid fourth-order derivative terms by intro- 
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ducing two auxiliary scalars. Then, the linearized curvature scalar becomes a scalaron, in- 
dicating that all linearized equations are second order. Using the positivity of gravitational 
potentials and S-deformed technique [22] allows us to prove the stability of f{R)- AdS black 
hole if the scalaron mass squared {m\) satisfies the Breitenlohner-Freedman (BF) bound 
{m\ > m| F = —^2) [23J. In order to confirm the stability of /(i?)-AdS black hole, we 
have to realize that a practical tool for testing stability in all those cases is a numerical 
investigation of quasinormal spectra k = k,R — ikj [13], [25]. We have checked that there 
is no exponentially growing mode (tachyon instability) in anti de Sitter space if the scalaron 
mass squared satisfies the Breitenlohner-Freedman bound (m^ > m^). For m\ < m| F , 
however, we have found a growing mode, showing the tachyonic instability of /(-R)-AdS 
black hole. 



2 Perturbation of f(R)-AdS black holes 



(2.1) 



Let us first consider f(R) gravity without any matter fields 



with k 2 = 8ttG. Since all relevant equations derived from (12 .ip are given by Ref. [21], we 
do not reproduce them here for saving the space. The point to remember is that there exist 
fourth-order linearized equations around the /(-R)-AdS black holes. In order to avoid it, we 
introduce two auxiliary fields and A which may allow to rewrite the action (12. ip as [IB] 

1 



S. 



2/t 2 



${R-A) + f{A) 



2.2 



Varying for the fields (f> and A lead to the following equations: 



R = A, <j> = f{A), 



(2.3) 



where ' denotes differentiation with respect to its argument. Note that imposing (12. 3p on 
the action (12. 2p recovers the original action (12. ip . On the other hand, the equation of 
motion for the metric tensor can be obtained by 



f'(A)R 



(^V 2 -V M V,)f(A)=0. 



(2.4) 
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In deriving the above equation, we used two relations in (12. 3p . Considering a constant 
curvature scalar R = R = A together with = f'(A) = const, Eq.f l2.4p becomes 



f'(A)R, u - \g^f{A) = 0. (2.5) 



Taking the trace of (12 .5p leads to 

n = nf s 4A " (2 ' 6) 

which determines the positive, negative and zero curvature scalar by choosing a specific 
form of f(A). 

At this stage, we wish to point out the difference between R = A in (I2.3P and (12. 6p . 
The former is the strong constraint, while the latter is satisfied only in the background of 
constant curvature scalar. In general, R is determined by the trace equation 

-3V 2 f(A) + 2f(A) 

R= — m — ■ < 27 > 

Substituting (12. 6p into (12. 5p . one finds the Ricci tensor which determines the maximally 
symmetric Einstein spaces including Minkowski space 

A /(i?)-AdS black hole solution is given by a spherically symmetric form 

ds 2 = g^dx^dx" = -e u{r) dt 2 + e~ u(r) dr 2 + r 2 (d9 2 + sin 2 6d^ 2 ) (2.9) 

with 



e"M = 1 - — - ^r 2 . (2.10) 

r 3 V ' 

We focus on asymptotically AdS4 spacetime with A^ = — J| < and consider the back- 
ground of 

R = A = 2f(A)/f'(A) < 0, (2.11) 

with f(A) < 0. 

Now we introduce the perturbation around the constant curvature black hole to study 
stability of the f(R)- AdS black hole 

9iiv = gp.v + h^ u . (2.12) 
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Hereafter the background quantities are denoted by the "overbar" notation. In this back- 
ground, we define the perturbations as 



R+6R{h) = A + SA, f(A)=f(A) + f'(A)8A, f'(A) = f'(A) + f"(A)5A. (2.13) 



The first relation implies the replacement of linearized curvature scalar by the linearized 
scalaron as 

8R(h)=R-R->8A = A-A. (2.14) 
Also the linearized equation to (12. 4p takes the form 

'f»(A)f(A)-r(Ay 



5R^(h) - A A h 



V' 2 {A) 



SA 



f"(A) 



(^V 2 -V M V 1/ )(JA = 0, 



(2.15) 



where the linearized Ricci tensor 5R^ v (h) is given by 
SR, u (h) 



1 

2 



(2.16) 



It is important to note that taking the trace of (I2.15P with g^ u together with (I2.14p leads 
to the linearized "scalaron" equation as 

(2.17) 



(2.18) 



V 2 - m 2 A )5A = 0, 
where the scalaron mass squared m 2 A is given by 



a f' 2 (A)-2f(A)f"(A) f'(A) 4 
mA ~ 3f(A)f'(A) ~3f"(A) 3 A ' 



which was already known as (97) of Ref. [2] in de Sitter spacetimes. Plugging ( I2.17P into 
(I2.15P and rearranging the terms, we obtain the linearized equation 



5R^(h)-A A h 



f"(A) 



f' 2 (A) + f(A)f"(A) 
Qf' 2 (A) 



5A. 



(2.19) 



Since the mass dimension of the scalaron is two ([SA] = 2), it would be better to write the 
canonically linearized equations by introducing a dimensionless scalaron 5 A defined by 
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(2.20) 



Here [f'(A)} = and [f"(A)\ = -2. 

Finally, we arrive at two linearized equations 

(V-m 2 W = 0, 



SR^(h) - A A V - V„V V 6A 



A, 



g^dA = 0, 



which should be solved to carry out the stability analysis of /(-R)-AdS black hole. 



(2.21) 
(2.22) 



3 Stability analysis of f(R)- AdS black hole 

The metric perturbations are classified depending on the transformation properties 
under parity, namely odd (axial) and even (polar). Using the Regge- Wheeler [25], and 
Zerilli gauge [27] , one obtains two distinct perturbations : odd and even perturbations. 
For odd parity, one has with two off-diagonal components ho and h\ 



h° 



( 











ho(r) \ 
h x (r) 




-ikt 



sin 9 



dpL 
dO 



(3.1) 



V h {r) hir) 0/ 



while for even parity, the metric tensor takes the form with four components Hq, Hi, H2, 
and K as 



( H (r)e u ^ Hi(r) 

#i(r) H 2 (r)e~ v ^ 

r 2 K{r) 

\ 









\ 



e~ lkt PL 



(3.2) 



sin 2 9K(r) J 



where pi is Legendre polynomial with angular momentum L and e u ^ was given by f)2.10p . 
In order to achieve the stability of f{R) black hole in asymptotically flat spacetimes, we 
have used the result for the stability analysis for the massive BD theory [19]. However, the 
present situation is quite different from f(R) black hole, because we are going to perform 
the stability analysis of /(-R)-AdS black hole. 
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3.1 Odd-perturbation 

For odd-parity perturbation, its linearized equation takes a simple form as 

SR^(h) - AaV = 0. 



(3.3) 



Using three equations of 8R t{p , SR rtp , and SRg^ in (13 .3p together with (I2.16P and (13. ip . one 
finds the following equation for h\\ 



d 2 h x 
dr 2 



, 2\ dh x 
r J dr 



v » + 2u' 2 + k 2 e 



„2-2v 



_„(L-l)(L + 2) 



/ii = 0. 



(3.4) 



Introducing the tortoise coordinate r* = J e u ^dr which maps r G [r + ,oo] into r* E 
[—oo,0] and a new field Q defined by 



r 

(13. 4 p leads to the Regge- Wheeler equation as 

ci 2 Q 



dr 



*2 



+ 



Q = 0. 



(3.5) 



(3.6) 



6m 

.3 



(3.7) 



Here the Regge- Wheeler potential is given by [26l [28] 

Tr r \ (■, 2m Ayi 2 \ \ L(L + 1) 
V RW {r) =1 —r ) 2 

It is well known from [28] that the effective potential can be positive definite outside the 
SAdS black hole when using S'-deformed technique [22J. This implies that the /(i?)-AdS 
black holes are stable against odd-perturbation because (13.61) is equivalent to that of the 
SAdS black hole in Einstein gravity. 



3.2 Even-perturbation 



We consider the linearized equation (I2.22p fully to perform even-parity perturbation. Con- 
sidering the metric perturbation (13. 2p and the scalar perturbation given by 



with Y^f(9,ip) spherical harmonics, the equation of either 

H 2 = Hq — 2—. 

r 



(3.8) 

or (ipip) component yields 

(3.9) 



From five remaining equations of (tr), (rr), (tt), (t9), and (r9) components, one obtains 
the constraint equation 

{Z 2 - 2e v + reV} H + ^2k 2 r 2 e~ v + 2e v + reV + ^r 2 eV 2 - Z 2 + 2A^r 2 1 A' 



- { 2ikr + ^}H X + |2L 2 - Ae u - 2k 2 r 2 e~ u + ^r 2 eV 2 - 2A A r 2 | ^ = (3.10) 

with L 2 = L(L + 1). In deriving this expression, we have used u" + (v 1 ) 2 + 2v' / r = —2k^e~ v . 
At this stage, we introduce the tortoise coordinate (r*) and a new field defined by 



where 



. . A(A + l)r 2 + 3Amr + 6m 2 , . . i(—\r 2 + 3Amr + 3m 2 ) 
<?( r ) = *n — rr~\ ' h (r) - 



r 2 (Ar + 3m) ( r - 2m - ^r 3 )(Ar + 3m) ' 

iv^ 1 

P(r) = — 3 , A = -(L-l)(L + 2). (3.12) 

r — 2m Y r ^ 

Manipulating (tr) and (tO) component equations in f)2.22p . we arrive at the Zerilli's equa- 
tiorfl 



d 2 M 
dr* 2 

where the Zerilli potential is given by [271 EH] 



V z 



M = 0, (3.13) 



/ 2m A A 2 

Vz(r) = 1 — r 

V r 3 



2A 2 (A + l)r 3 + 6X 2 mr 2 + 18Am 2 r + 18m 3 - 6A A m 2 r 3 



r 3 (Ar + 3m) 2 

(3.14) 

The Zerilli potential Vz is always positive for whole range of — oo < r* < (see Fig. 9 
of [22]), which implies that the even-perturbation is stable, even though the scalaron ip(5A) 
is coupled to making the even-perturbation. 



1 At this stage, one may ask the question of "can one find the single metric components (e.g. H\ or ip) 
from the master variable M?" . The answer is obviously "No" . Actually, one has /i M „ with 10 components 
in Einstein gravity. A masslcss graviton has 2 degrees of freedom (DOF). Upon the Rcggc- Wheeler gauge- 
fixing, one has 6. One DOF is given by the odd-perturbation combined by 2(hi and /12) and the other is 
the even-perturbation combined by 4 (Ho, Hi, H2, K). In f(R) gravities, one starts with 11 (10 from 
and 1 from ip), indicating 3 DOF (masslcss graviton and scalaron) totally. Hence, the scalaron should be 
coupled to giving the master variable M in the even-perturbation. 
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3.3 Scalar-perturbation 



The main difference between Einstein and f(R) gravities is the appearance of scalaron. 
Hence, the stability analysis of scalaron equation will be considered to be the main part of 
the present work. In order to obtain the second-order differential equation for the scalaron, 
we first recall that the scalaron perturbation is given by (13. 8p . In this case, Eq. (12.2 lj ) 
becomes 

L(L + 1) 



eV 



r r 



m\ — e u k 2 



+ 



eV 



+ 



ip = 



(3.15) 



Using the tortoise coordinate (r*), the above equation reduces to the Schrodinger-type 
equation 



dr 



*2 



k 2 -V A 



i' = o, 



where the scalaron potential Va is given by 



V A (r) 



2m 

r 



2 

— r 
3 



L{L + 1) 2m 



m. 



(3.16) 



(3.17) 



r* r J 3 

Here, the second term in square bracket is the usual scalar term with spin zero [in general, 
—2m(s 2 — 1) for s spin-weight of a perturbing field], and the third arises from AdS4 asymp- 
tote. Importantly, the last term in bracket shows the feature of a massive scalaron arisen 
from f(R) gravity. This term could not be eliminated using the S'-deformed technique 
because it is not an element of the metric function e vtj "\ Furthermore, it is well known that 
there exists the BF bound [23] for the Klein-Gordon type equation (I2.17P 



which implies an inequality 



m\ > m 2 BF = - A 



f(A) 25 
3f"(A) ~ U 2 ' 



On the other hand, we rewrite the scalaron potential (13.171) 



V A (r) 



2m 
r 



A 



A 2 

r 



L(L 



1) 2m 
— + —r 



6 f(A) 
i 2 3f" {A) 



(3.18) 



(3.19) 



(3.20) 



which shows the asymptotic positivity (V^(r — > oo) — > 0) when the following condition is 
achieved 

f'(A) 6 

(3.21) 



6 

3f"{A) ~ "i 2 ' 
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Hence, we allow the negative mass squared in A0IS4 spacetimes which belongs to the region 

25 f'(A) 6 9 „ 2 , 

< - v L < y < ml < . (3.22) 

U 2 ~ 3f"(A) ~ P AP ~ A ~ P K ' 

We note that the left bound comes from the BF bound, while the right bound arises from 
the asymptotic positivity condition of Va{t —> 00) —> 0. Accordingly, we divide m\ into 
four regional m\t 2 < m| F £ 2 = -9/4, - 9/4 < m\t 2 < -2, - 2 < m\t 2 < 0, and 
m 2 A i > 0. The boundaries are m 2 A l 2 = 0, —2, —9/4. As is shown in Fig. 1, there are four 
types of graphs including three boundaries, depending on the mass squared m\. This may 
imply that there is no ghost instability even for the negative mass squared is allowed, if m\ 
satisfies the bound (13.221) . This is surely a new feature of stability condition for /(i?)-AdS 
black holes. On the other hand, we expect that for m\ < m| F , there will be the tachyonic 
instability. 



4 Quasinormal modes 

It is well known that a practical tool for testing stability of all kinds of black holes is 
a numerical investigation of the quasinormal spectra [T3|. Hence, we wish to confirm the 
stability of /(i?)-AdS black holes by computing quasinormal modes. Considering ty(t, r*) = 
e~ lkt ty(r*), the boundary conditions for quasinormal modes propagating in the /(-R)-AdS 
black hole background are given by 

{ e -ikr« for r* -»• — 00, 
(4.1) 
0, for r* -)> 0, 

where ^(r*) represents Q (odd), M (even), and ip (scalaron) in (13. 6j) , (13 . 1 3j) . and (I3.16p . 
respectively. The frequencies of the quasinormal modes are given by complex numbers, 



2 The fourth region m\£ 2 > corresponds to "tachyon-free condition" required by the stability of de 
Sitter space of (/' 2 - 2//")/3/7" > E H3 130] and Minkowski space of f'/3f" > [18,. This implies 
that /' > and /" > for stability. However, the stable region (tachyon-free case) could be extended to 
negative mass squared in the AdS space: m\£ 2 > — 9/4[ (|3.19[) ]. which implies that f"(A) < 0, for f'(A) > 
in AdS space. This is in contrast with other perturbation analysis: the Dolgov-Kawasaki instability with 
/" < in cosmological perturbations [31] , graviton and scalar propagations in the Minkowski [18l [32] • 
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2 = 1: for m 2 A i 2 = case, 



Figure 1: Scalaron potential Va graphs as function of r with fixed » 
it shows the massless scalaron in AdS spacetimes, while for m 2 A £ 2 = —2 case, it approaches 
zero as r — > oo (asymptotic positivity). For two cases of m 2 A £ 2 = —9/4 and —5/2, their 
potentials approach negative infinity as r — > oo. 

k = kR — iki, and the modes with kj > (kj < 0) corresponding to decaying (growing) 
ones, since the time dependence is given by e~ lkt = e~ 
least one growing mode, the black hole is unstable. 
Now, (13. 6p . (I3.13p . and (13 . 1 6[) could be written as 



~kjt ~—ihf>t 



e R . We note that if there is at 



drl 



(k - V)V(r t ) = 



(4.2) 



with V denoting Vrw, Vz, and Va- Multiplying this equation by a complex conjugated 
function ^(r*) and integrating it, we have 



-4k l^(-oo) 



d 2 V(r, 
dr 2 



dr* 



{k 2 - V) |#(r*) 



dV(r* 



dr* 



(k 2 -V) |tf(r, 



(4.3) 



where we used the boundary condition (14.11) . The imaginary part of (14. 3 j) takes the form 



Im(7) 



l^(-oo) 



2k! 



dr* \^f(r* 



0. 



(4.4) 



which implies that either k^ = or kj > 0, so that the growing modes cannot oscillate as 
was pointed out in [251 US]- That is, the unstable mode is defined by the condition of 



unstable mode k^ = 0, ki < 0. 



(4.5) 
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in the quasinormal-mode approach. Next, the real part of (14.31) is given by 



Re(J) =kj |^(-oo) 



dr* 



d^(r 



0. 



(4.6) 



which means that the unstable mode does not exist if the potential V is positive definite. We 
confirm that there is no unstable modes for the even-perturbation since Vz is always positive. 
For the odd-perturbation, Vrw becomes negative near the horizon, but the stability may 
be achieved by the S-deformation technique. On the other hand, it turned out that there 
is no growing modes in odd- and even-quasinormal mode computations [2"8"| I3"3"]. Hence, we 
do not need to compute quasinormal frequencies for the metric perturbations again. 

However, we remind the reader the instability of the scalar perturbation depends on 
the scalaron mass squared, m\. Requiring Va to be positive corresponds to the bound of 
m 2 A > —2£~ 2 . Thus, one expects that there is no exponentially growing mode for this case. 
Also, if m A satisfies the BF bound (m| F < m 2 A < —2£~ 2 ), there is no tachyon instability. 
Finally, if m 2 A < m| F , there will be tachyonic instability. In order to show it explicitly, we 
calculate quasinormal frequencies of the scalaron numerically. 

To this end, Eq. (14. 2p should be rewritten in the ingoing Eddington coordinates with 
v = t + r*, following the original work of Ref.JM] for large Schwarzschild-AdS black holes: 



,d 2 V 
dr 2 



+ {v'e v - 2ik) 



d^ 

dr 



= 0, 



(4.7) 



where ^(t, r*) 



-ikv 



and V = e v V . Then, changing coordinate variable r into x 



the equation (14. 7p becomes 



. .d 2 ^ 
S{x) lx^ + 



t(x) d^ 



u{x) 



X 



x dx {x *^+) 



= 0, 



where 



s(x) 



x 4 e u 



x j 



X 



— [x\ + x + x + (1 + x 2 + )x 2 ] = ^ s n (x - x + ) n , 

X j 

3 



n=0 



t(x) = - 2ikx 2 - 2x 3 + —(1 + x 2 + )x 4 = ^t n 



X 



X. 



n=0 



u(x) ={x — x + )V = u n (x — 



X A 



(4.8) 

(4.9) 
(4.10) 
(4.11) 



n=0 
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Table 1: Frequencies of quasinormal modes of large /(-R)-AdS black holes for the scalaron 





10 


50 


100 


k R 


~ki 


k R 


h 


k R 


kj 





18.6070 


26.6418 


92.4937 


133.193 


184.953 


266.386 


-1/4 


18.0708 


25.6542 


89.8263 


128.258 


179.619 


256.515 


-1/2 


17.500 


24.6056 


86.988 


123.018 


173.945 


246.036 


-3/4 


16.888 


23.4827 


83.941 


117.407 


167.85 


234.813 


-1 


16.222 


22.2661 


80.629 


111.327 


161.23 


222.653 


-5/4 


15.486 


20.9261 


76.968 


104.630 


153.91 


209.260 


-3/2 


14.64 


19.42 


72.81 


97.067 


145.6 


194.13 


-7/4 


13.66 


17.63 


67.88 


88.15 


135.7 


176.3 


-2 


12.3592 


15.3307 


61.4106 


76.6646 


122.797 


153.330 


-9/4 


9.9 


12 


49 


58 


98 


116 


-5/2 





-8.64855 





-32.5878 





-59.5519 


-11/4 





-52.5789 





-225.272 





-467.935 


-3 





-126.996 





-536.309 





-1035.61 



Here x + = £/r + and k = k£, and the boundary condition (14. ip takes the form of ^(x — > 
0) —¥ 0. Considering a series solution 

oo 

§ = ^a n (x-a; + ) n , (4.12) 

n=0 

the boundary condition at infinity yields an algebraic equation 

oo 

^a n (-x+r = 0. (4.13) 

n=0 

Solving this, one could find (discrete) quasinormal modes. In Table [TJ we list quasinormal 
modes for the scalaron (L = 0) obtained by truncating the series after 140 terms. The first 
row of m 2 A £ 2 = indicates the quasinormal mode for the massless scalar, which is consistent 
with Ref. [31]. The last three rows in Table 1 shows clearly that there exist unstable modes 
of kj < with k R = for the scalaron mass (m 2 A £ 2 = —5/2 < —9/4) below the BF bound. 
This corresponds to the tachyonic instability. For —9/4 < m 2 A £ 2 < —2, there is no unstable 
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mode, showing that there is no tachyon instability. This confirms the presence of the BF 
bound even for f(R)- AdS black holes. 

5 Discussions 

We have investigated the stability of /(i?)-AdS (Schwarzschild-AdS) black hole obtained 
from the f{R) gravity. Even though the Schwarzschild-AdS black hole is known to be 
stable in Einstein gravity, checking the stability of /(i?)-AdS black hole is a nontrivial task 
because the linearized Einstein equation is fourth order in f(R) gravity. In order to resolve 
this difficulty, we have translated the fourth-order equation into the second-order equation 
by introducing auxiliary scalar fields A and 0. In this case, the linearized curvature scalar 
5R(h) becomes a massive scalaron 5A, showing that all linearized equations are second 
order. 

The stability on the metric perturbations remains unchanged, confirming that the odd 
(even) perturbations lead to the Schrodinger-type equation with the Regge- Wheeler (Zerilli) 
potentials in asymptotically AdS 4 spacetimes. Actually, this corresponds to the case in 
Einstein gravity, even though the even-perturbation contains the scalaron in addition to 
Ho, H\, H2, and K. For the odd-perturbation, we have needed the S'-deformation method 
to prove the stability of /(i?)-AdS black hole because the corresponding potential becomes 
negative near the event horizon. 

The main difference comes from the linearized scalaron equation because the scalaron 
plays the role of a massive scalar, which is physically propagating on the /(i?)-AdS black 
hole background. We have shown that /(i?)-AdS black hole is stable against the scalaron- 
perturbation if the scalaron mass squared m 2 A satisfies the BF bound. Especially for —9/4 < 
m\ < —2, there is no unstable mode (exponentially growing mode), showing that there is 
no tachyon instability. This confirms the presence of the BF bound even for /(i?)-AdS black 
holes. This is consistent with graviton and scalar propagations AdS [35] spacetimes, but 
contrasts to other perturbation analysis: the Dolgov-Kawasaki instability with f"(R) < in 
cosmological perturbations [21], graviton and scalar propagations in the Minkowski [TBI I3"2"] . 
and de Sitter spacetimes [29], [30] . 

Finally, we have found the tachyonic instability for m\ < m| F . This is confirmed by 
computing quasinormal frequencies of the scalaron for large /(-R)-AdS black hole. The last 
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three rows in Table 1 shows clearly that there exist unstable modes of kj < with A;^ = 
for the scalaron mass (m\£ 2 = —5/2 < —9/4) below the BF bound. This indicates the 
tachyonic instability. 
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